Excited state TBA and functional relations in spinless Fermion model 
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The excited state thermodynamic Bethe ansatz (TBA) equations for the spinless Fermion model 
are presented by the quantum transfer matrix (QTM) approach. We introduce a more general 
family called T-functions and explore functional relations among them (T-system) and their certain 
combinations (Y-system). From their analytical property, we derive a closed set of non- linear integral 
equations which characterize the correlation length of (cjci) at any finite temperatures. Solving these 
equations numerically, we explicitly determine the correlation length, which coincides with earlier 
results with high accuracy. 



Thermodynamics of ID quantum integrable systems 
have been discussed by the thermodynamic Bethq-ansatz 
(TBA) methodEl based on the string hypothesis.ETtl Ther- 
mal quantities are determined by a set of non-linear inte- 
gral equations called TBA equations. As an alternative 
method, the quantum, transfer matrix (QTM) approach 
has been proposedJa~c3 It utilizes the Trotter formula and 
reduces the calculation of the free energy to the eigen- 
value problem of the QTM on a fictitious system. A set 
of auxiliary functions, including the QTM itself, satisfy 
functional relations. Under suitable choice, these func- 
tions have a property called ANZC (Analytic, NonZero 
and Constant asymptotics) in certain strips. This admits 
a transformation of the functional relations to non- linear 
integral equations (NLIE) which characterize the free en- 
ergy. 

Though each method reproduces the correct free en- 
ergy, there seems to be no relation between the QTM 
and the TBA. Recently, however, we recognize the NLIE 
are identical to the. TB.A equations by utilizing the fu- 
sion hierarchyOO'LiltS The members of fusion hierar- 
chy (T-functions) and the combinations of them (Y- 
functions) satisfy functional relations called T-system 
and Y-system, respectively. (See ref. |27],^8| for general re- 
lations between T and Y-systems.) By selecting these 
functions such that they are ANZC in appropriate strips, 
we derive the NLIE which are identical to the TBA equa- 
tions. Moreover considering the sub-leading eigenvalues, 
we derive systematically the "excited state TBA" equa- 
tions, which is difficult within the string hypothesis. 

We apply this procedure to the spinless Fermion model. 
The physical quantities such as the free energy and the 
correlation length have been already derived by more 
practical choice of the auxiliary functions in ref. Eq. Our 
present interest, however, is the derivation of the TBA 
and excited state TBA equations by utilizing the fusion 
hierarchy and to confirm the consistency between the re- 
sultant equations and the earlier ones. The Hamiltonian 
on the periodic lattice of size L is 
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where we consider the model in the repulsive critical re- 
gion 0< A = cos6>< 1, < t and introduce the parameter 
Po as po = ir/9(po > 2). The commuting QTM on a 
fictitious system of size TV (Trotter number, N S 2Z) has 



been defined by. the Fermionic i?-operator and its super- 
transposition.c3 The eigenvalue is expressed as 



T 1 (u,v) = <f,+ (v)<f>-(v-2i) 
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Here e = (-l) N / 2 + N * and N c = {0, 1, • • • , N/2} is the 
quantum number counting the holes on odd sites and 
particles on even sites, {wj} is a set of solutions to the 
Bethe ansatz equation (BAE): 



</>+(v)<l>-(v - 2i) 
<j>-{v)cj> + {v + 2i) 
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Q(wj - 2i) 
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The largest eigenvalue T^\u,v) and the second largest 
eigenvalue t[ 2) (u,v) lie in N c = n/2 and iV c =N/2-l, 
respectively. Hereafter we often omit the u variable in 
Ti(«, v). The free energy per site / is represented by 

(3) 

The correlation length £ of (cjCi) can be also given by 
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the ratio of T^\v) to t[ 2) {v) 
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To evaluate T±(v), we consider a more general family 
called T-functions: 

n 

T„_i(u, v)= J2 V- (« - i(n + 2 - 2j)) 

3 = 1 

4>+{v — i(n — 2j))Q{v + in)Q(v — in) 



Q(v + i(2j - n))Q{v + i(2j - n - 2)) 
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where we set T-\{x) = 0. 

For any v £ C and integers n > m > 1, the T-functions 
satisfy the following relations, 
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l T m _i(u + in)T m _i(w - in). (6) 



From now on we consider the case po S Z>3. In this case 
T-functions satisfy further functional relation, 

T P0 (v) = sT po _ 2 (v) + (-1)^(1 + e Po )To(v + ip ). (7) 

We call the above two functional relations the T-system. 
The proof of them are direct by using (^) and the peri- 
odicity of the T-functions 

T n _i(u) = T n - 1 {v + 2 Po i). (8) 

Let us define a set of combinations among T-functions as 



e 3 + Yj[v) 
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K(v) = (-1) 



T 



Po-2 



T (u + ip ) 
where 1 < j < p — 2, 1q( v ) 



and y-i(w) = oo. We 



find that {^(v)}^! and if(u) satisfy following finite 
set of functional relations, which we call the Y"-system. 

- <*i- x +Y i ,(v))(E j+1 +Y, 



Yj(v + i)Y s (v - i) = (e^ 1 + 

1 1 v = (e + X(v))(e w 
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+ + Y Pa - 1 (v) = {e + K{v)){e P0 +sK(v)), (10) 
Jif(« + - i) - e P0 + Y P0 _ 2 (v), 

This can be proved by the T-system (|TJ) and (Q). Let 

T^^v) denote T„_i(w) constructed by the BAE roots 
which is relevant to the fc-th largest eigenvalue and 
Yj k \v), K^(v) be the F-functions constructed from 

{^n-i} as m ©• Using the property </>±(v ± 2poi) = 
(-l^ViO) and Q(« + 2p i) = (-l) w,/a Q(«), we can 
easily show for po € 2Z 



7f^) = o, y p ( 2 l 2 w-o. 
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To proceed further, one clarifies the analyticity of the 
Y"-functions by numerical studies on the T-functions, 
keeping the Trotter number N finite. From them, af- 
ter suitable modification we confirm the V-functions 
have the ANZC property in the strip Imw S [—1,1]. 
We call this "physical strip" . Then we can transform 
the F-system to a closed set of NLIE in the following 
way. First we take the logarithmic derivative and per- 
form the Fourier transformation on both side of the first 
and third equation in (|To|) . Second by Cauchy's the- 
orem, the Fourier mode for the logarithmic derivative 
of Yj(v) (K(v)) is expressed by those of + Yj+i(v) 
and e^ 1 + Y}_i(u) (e Po + Yp n _2(v))- Finally performing 
the inverse Fourier transformation and integrating over 
V, we derive the desired NLIE. In these NLIE, the Trotter 
limit N — > oo can be taken analytically. The eigenvalue 
T^ fe ' (v) are determined after performing such transfor- 
mation on 



Ti(« + i)Ti(> - i) = (e + Y^v^Toiv + 2i)T {v - 2i), 



For the largest eigenvalue in the sector iV c = N/2, all 
functions and functional relations are equivalent to those 
of the XX Z model in the critical regime.Eil Therefore one 
concludes that the resultant NLIE which characterize the 
free energy for the spinless fermion model .are identical 
to the TBA equations for the XX Z models-El 

Let us consider the second largest eigenvalue in the 
sector N c — N/ 2 — 1 . As is mentioned in ref . we find 
that two pure imaginary eigenvalues which are complex 
conjugate each other, are degenerate in magnitude. All 
the BAE roots for these two eigenvalues are real. The 
two distributions of the BAE roots are symmetrical with 
respect to the imaginary axis. Hereafter we consider the 



case ImTf^O) > (the case ImT^^O) < also can be 
discussed in similar way). According to the numerical 
studies, we have the following conjecture which is quite 
different from that in XX Z model£il 



(2), 



(2) 

Conjecture 1 T^l^v) has one real zero Cn— l f or { n G 
2Z|2 < n < po — 1} and two real zeros Cn-i an d ~ Cn-i 
for {nE2Z + l|2<n< po/2}. All other zeros are out 
of the physical strip. 

Here Cn-ij Cn-i > 0- For example we depict the location 
of zeros in Fig. [l] for po = 10. 
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FIG. 1. Location of zeros for T^2i(u,v) for u = —0.05, 
po = 6, N = 12. There exists one real zero for n = 2 and 
n = 4, while for n — 3 two zeros are at ±oo. Zeros on the 
imaginary axis for n = 5 are sextuple roots. 

Y^ and K^> have the asymptotic values 
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for j € 21 
for j e 2Z - 
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for po € 2Z 
itan# (v — » 00) for po S 2Z ■ 
— itan0 (y — > — oo) for po S 2Z ■ 



(13a) 



(13b) 



From conjecture 1 and above asymptotics, the factors of 

-iy +1 + Y?\(v) are ANZC in 



(1^) the rhs in (10) such as 
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Imv £ [— e,e] (e <C 1)- For the lhs, we need to modify 
Y} 2 \v) ({3 G 2Z|2 < j < p - 3} or {j e 2Z + 1|1 < 
j < Po/2}) and i^ 2 )(i;) (p G 2Z + 1) such that they are 
ANZC in the physical strip. For example, we take the 
case po — 6. The following modification, 
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is sufficient. Note that the factor e ™th(™/4)/ Po ig in _ 
eluded to compensate the singularity caused by Y^(v) 

and y£ 2 \v) tending to zero as e~^ v ^ Po at v = ±00. Let 
us rewrite the ^-functions in the Trotter limit as 



m(v)= lim y} 2 \u n ,v), k(v)= lim K {2) (u N , v). (f 4) 



After taking the logarithmic derivative and performing 
the Fourier transformation, we derive the NLIE obeyed 
by rjj and k. 



In 771 (u) = 
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Si * ln((l + r/ 2 )/ii)(u) 



7ri for po 5: 5 

Tri - f uth^ for Pa = 6 

ln{th|(«-C 2 )thf(u + C2)} forp >7 



ln % («) = si * ln(l - r^iXl - + ln{th-(« - 0_i)th-(u - 0+0} + *i for {j G 2Z|2 < j < p - 3}, 

lnrj^u) = si * ln((l + ?7j_i)(l + rjj + i)hj)(v) 
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V Po -2{v) = 0, k(v) = 1 

lmj Po - 2 {v) =si* ln(l + ry Po _ 3 )(l - k 2 )(") 
L ln/c(t>) =si* ln(l - r) Po - 2 )(v) + ln{th|(« - C Po - 
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where 



hj(v) 



exp 
1 
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fbr{j G 2Z- 
othcrwisc. 



The symbol * denotes the convolution 

f*g(x)=l f{x-y)g(y)dy, 



Po-4 po -j 
2 ' 2 J 



and si(u) = l/4chSu. Here the integration cons tants 



have been fixed from the asymptotic values ( 13a ) and 

( |13b| ). In addition to above equations, we need to im- 

t c ?\ 

pose the consistency condition coming 



and ri 2) (-Cj 



0. From 



and fll4!), this leads to 
T]j(Q ± i) = 1 for {j G 2Z + 1 1 1 < j < po - 2} and 
ife(0±i) = ife(Cj±i) = -1 for {j G 2Z|1 < j < Po/2-1}. 
Explicitly they read 
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< i < f } 



po j 



for {j G 2Z + l|j - , 
for {je2Z+l|f <j<p -3} 

for po € 2Z + 1 and p 7^ 3, 
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where the convolutions should be interpreted as 
si *g(C + i) = p.v 



9{x) d X \ + Uc). 



-co 4ishf (C-x 



Here p.v. means the principal value. Then the Trotter 
limit of T-J (v) can be expressed from ( |l2] ) as 

lim r 1 (2) (« Ar ,0)=]nth5Ci 

N^oo 4 



cfosi(i>) ln(l — 



+t/3 / dusi(w) 



sin i n . it . 

a cos a H — 1. 

ch9v-cos6 T 2 



Especially for po — 4, the zero £i and Y-functions 771 (i>) 
are explicitly determined by using the fact (11), which is 
a novel feature in the present approach. cB 

Finally we obtain the correlation length (Q) as 



- = - lnth ^Ci - f° dv 8l (v) In ( ^— ] , (17) 

where the function 77^ (u) characterizing the free energy 
has been calculated in ref. E 21. We evaluate (|l5j) and 
( |l6| ) numerically and explicitly determine the correlation 
length in Fig. 0. 
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FIG. 2. Ratio of the correlation length and inverse temper- 
ature for po — 3, 4, 5, 6 and 7. 

These results are consistent with those in ref. Espe- 
cially in the low temperature limit they agree with known 
expression from CFT 



lim &(/3)//3 = 



t sin ( 
26 



4(tt - 6) 
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